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Abstract 



Without additional resources, it is often impossible to transform one entangled 
quantum state into another with local quantum operations and classical communica- 
tion. Jonathan and Plenio [Phys. Rev. Lett. 83, 3566(1999)] presented an interesting 
C<") , example showing that the presence of another state, called a catalyst, enables such a 

transformation without changing the catalyst. They also pointed out that in general 
it is very hard to find an analytical condition under which a catalyst exists. In this 
paper we study the existence of catalysts for two incomparable quantum states. For 
, the simplest case of 2 x 2 catalysts for transformations from one 4x4 state to another, 

a necessary and sufficient condition for existence is found. For the general case, we 
give an efficient polynomial time algorithm to decide whether a k x k catalyst exists 
O ^ for two n x n incomparable states, where k is treated as a constant. 

+!> ■ 

Index Terms — Quantum information, entanglement states, entanglement trans- 
formation, entanglement catalysts. 

> 

X' 

Entanglement is a fundamental quantum mechanical resource that can be shared among 
spatially separated parties. The possibility of having entanglement is a distinguishing fea- 
ture of quantum mechanics that does not exist in classical mechanics. It plays a central role 
in some striking applications of quantum computation and quantum information such as 
quantum teleportation pQ, quantum superdense coding |2j and quantum cryptography 
As a result, entanglement has been recognized as a useful physical resource (!]. However, 
many fundamental problems concerning quantum entanglement are still unsolved. An 
important such problem concerns the existence of entanglement transformation. Suppose 
that Alice and Bob each have one part of a bi-partite state. The question then is what 

'This work was partly supported by the National Foundation of Natural Sciences of China (Grant Nos: 
60223004,60496321,60321002, and 60305005). 
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other states can they transform the entangled state into? Since an entangled state is 
separated spatially, it is natural to require that Alice and Bob can only make use of local 
operations and classical communication (LOCC). Significant progress in the study of en- 
tanglement was made by Bennett, Bernstein, Popescu and Schumacher p] in 1996. They 
proposed an entanglement concentration protocol which solved the entanglement trans- 
formation problem in the asymptotic case. In 1999, Nielsen j^| made another important 
advance. Suppose there is a bi-partite state \ipi) = Y17=i v^il*)^!*) B shared between Alice 
and Bob, with ordered Schmidt coefficients (OSCs for short) a.\ > ai > ■ ■ ■ > a n > 0, and 
they want to transform \tpi} into another bi-partite state \ip2) = Y27=i VAK)a|^)b with 
OSCs Pi > 02 > • • • > Pn > 0. It was proved that \ipi) — > \ip2) is possible under LOCC if 
and only if -< A</, 2 , where and X^ 2 are the vectors of ordered Schmidt coefficients, 
i.e. X l f >1 = (ai, . . . ,a n ), A^ 2 = (Pi, . . . ,P n ), ~< denotes the majorization relation [3IH|, i-e. 
for 1 < I < n, 



with equality when I = n. This fundamental contribution by Nielsen provides us with an 
extremely useful mathematical tool for studying entanglement transformation. A simple 
but significant fact implied by Nielsen's theorem is that there exist incomparable states \ipi) 
and \ip2) with both transformations \ipi) — > |^) an d \1jj2) — * \ipi) impossible. Shortly after 
Nielsen's work, a quite surprising phenomenon of entanglement, namely, entanglement 
catalysis, was discovered by Jonathan and Plenio [§]. They gave an example showing 
that one may use another entangled state |c), known as a catalyst, to make an impossible 
transformation — > possible. Furthermore, the transformation is in fact one of 
\ip) <8> |c) — * \4>) <g) \c), so that the catalyst |c) is not modified in the process. 

Entanglement catalysis is another useful protocol that quantum mechanics provides. 
Therefore to exploit the full power of quantum information processing, we first have to 
solve the following basic problem: given a pair of incomparable states \ipi) and \ip2) with 
IV'i) ~h 1^2) and 1^2) ~h l^i)) determine whether there exists a catalyst |c) such that 
IV'i) ® l c ) ¥ 1^2) 8> |c). According to Nielsen's theorem, solving the problem requires 
determining whether there is a state \c) for which the majorization relation A^ l8)c -< \ip 2 ® c 
holds. As pointed out by Jonathan and Plenio [H], it is very difficult to find an analytical 
and both necessary and sufficient condition for the existence of a catalyst. The difficulty 
is mainly due to lack of suitable mathematical tools to deal with majorization of tensor 
product states, and especially the flexible ordering of the OSCs of tensor products. In 
[0], Jonathan and Plenio only gave some simple necessary conditions for the existence 
of catalysts, but no sufficient condition was found. Those necessary conditions enabled 
them to show that entanglement catalysis can happen in the transformation between 
two n x n states with n > 4. One of the main aims of the present paper is to give 
a necessary and sufficient condition for entanglement catalysis in the simplest case of 
entanglement transformation between 4x4 states with a 2 x 2 catalyst. For general case, 
the fact that an analytical condition under which incomparable states are catalyzable 
is not easy to find leads us naturally to an alternative approach; that is, to seek some 
efficient algorithm to decide catalyzability of entanglement transformation. Indeed, an 
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algorithm to decide the existence of catalysts was already presented by Bandyopadhyay 
and Roychowdhury |l()j . Unfortunately, for two nx n incomparable states, to determine 
whether there exists a k x k catalyst for them, their algorithm runs in exponential time 
with complexity 0([(nk)l] 2 ), and so it is intractable in practice. The intractability of 
Bandyopadhyay and Roychowdhury 's algorithm stimulated us to find a more efficient 
algorithm for the same purpose, and this is exactly the second aim of the present paper. 

This paper is organized as follows. In the second section we deal with entanglement 
catalysis in the simplest case of n = 4 and k = 2. A necessary and sufficient condition 
under which a 2 x 2 catalyst exists for an entanglement transformation between 4x4 states 
is presented. This condition is analytically expressed in terms of the OSCs of the states 
involved in the transformation, and thus it is easily checkable. Also, some interesting 
examples are given to illustrate the use of this condition. The third section considers 
the general case. We propose a polynomial time algorithm to decide the existence of 
catalysts. Suppose l^i) and l^) are two given n x n incomparable states, and k is any 
fixed natural number. With the aid of our algorithm, one can quickly find all k x k 
catalysts for the transformation \ipi) — > ^2) using only 0(n 2fc+3 ' 5 ) time. Comparing to 
the time complexity 0([(nk)\] 2 ) of the algorithm given in JO]; for constant k, our algorithm 
improves the complexity from superexponential to polynomial. We make conclusions in 
section 4, and some open problem are also discussed. 

To simplify the presentation, in the rest of the paper, we identify the state \tp) = 
Yji=i w ffh the vector of its Schmidt coefficients (71,72, • • • ,7n)> the meaning will 

be clear from the context. 

2 A necessary and sufficient condition of entanglement catal- 
ysis in the simplest case (n = 4, k = 2) 

Jonathan and Plenio [9] has shown that entanglement catalysis only occurs in transfor- 
mations between n x n states with n > 4. In this section, we consider the simplest 
case that a transformation from one 4x4 state to another possesses a 2 x 2 catalyst. 
Assume \tpi) = cti, «3> CK4) and ^2) = (Pi, P2, P3, Pi) are two 4x4 states, where 



cui > ol 2 > a 3 > a A > 0, Yn=i a i = 1 , Pi > P2 > P3 > Pi > 0, and £i=iA = 1. 
The potential catalyst is supposed to be a 2 x 2 state, denoted by \<f>) = (c, 1 — c), where 
c E [0.5,1]. 

It was proved in jHj that if ~h but IV'i) <8> \4>) — ► ^2) &> \4>) then 



Oil < Pi, Oil + «2 > Pi + P2 , Oil + Ol2 + «3 < Pi + P2 + Pz 



(1) 



or equivalently, 



Oi2 + 013 + a 4 > Pi + p 3 + Pi, a 3 + a 4 < p 3 + Pi, a 4 > Pi. 



(2) 



Note that {a{\ and {Pi} are arranged in decreasing order, so we have 



Pi > 011 > "2 > P2 > P3 > 013 > a 4 > Pi 



(3) 
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These inequalities are merely necessary conditions for the existence of catalyst \<f>), and it 
is easy to see that they are not sufficient. In the following theorem we give a condition 
which is both necessary and sufficient. 

Theorem 2.1 There exists a catalysts \<f>) for two states (|^i), l?/^)) with \ip\) \if} 2 ), if 
and only if 

(01 + 02-/31, a 4 -0 4 \ . f (3i 0i -a ± 4 \ 

max — — — - ,1-- <mm — - — , -=-,1 — > (4) 

I 02 + 03 P3 ~ «3 J I ol\ + a 2 a 2 - 2 a 3 + a 4 J 

and (QJ) /10W. In addition, for any c 6 [0.5, 1] smc/i i/iai 

' ai + a 2 - A 1 «4 - /?4 1 . ^ . f A Pi -ax 0. 
max < , 1 — > < c < mm < , — , 1 



02+03 03 — "3 J I «1 + «2 «2 —02 «3 + «4 

|0) = (c, 1 — c) is a catalyst for (\ipi), \ip 2 })- 

Proof: Assume \ipi) ^2) but (8) \</>) -> |^ 2 ) ® |</>) under LOCC. From Eq. (8) in [§] 
we know Eq. (0) holds. So Eq. (j2j) and Eq. © hold too. 

A routine calculation shows that the Schmidt coefficients of \^i)\4>) and \ip 2 )\4>) are 

A = {a\c, a 2 c, a 3 c, a^c; q?i(1 — c), 02(1 — c),a 3 (l — c), 04(1 — c)} 

and 

B = {0 1C , 2 c, 3 c, 04c; 0i(l - c),0 2 {l - c),0 3 {l - c),0 4 (l - c)}, 

respectively. Sort the elements in A and B in decreasing order and denote the resulted 
sequences by > > ■■■ > and b^ 1 ' > > ■■■ > b^ . It is clear that aW = aic, 
a ( § ) = a 4 (l - c), &W = lC , and &( 8 ) = 4 (l - c). Since \ipi) \<f>) -» |V> 2 > ® Nielsen's 
theorem tells us that 

< b(i) (VI < ^ < 8) 

i=l i=l 

Since {0i} is ordered, and c > 0.5, thus 

0ic > 2 c > 3 c > 04C, 0i(l-c)>0 2 (l-c)>0 3 (l-c)>0±(l-c), 0iC>0i(l-c) (5) 

Now we are going to demonstrate that 

0ic > 0i(l -c)> 2 c > 3 c > 2 {1 - c) > 3 (1 -c)> /? 4 c > p A (l - c). (6) 

and consequently fix the ordering of B. The key idea is: the sum of the biggest / numbers 
in a set is greater than or equal the sum of any I numbers in this set. 

First, by definition of {o®} we have + > a\c + a 2 c. So Nielsen's theorem 
leads to &W + fr( 2 ) > a (X) + > a\c + 02c. From inequality (0), «i + "2 > 01 + 02, so 

_|- 5(2) > 0xc + 2 c, i.e. b^ > 2 c. Combining this with inequality ©, we see that the 
only case is b^ = 0i(l - c), &( 3 ) = 2 c and 0i(l - c) > 2 c. 
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Similarly, we have 

a (1 ) + + a® + o (4 ^ > a\c + a 2 c + a\{l — c) + 02(1 — c) = ai + a 2 . 
So it holds that 

5 (D + h m + 6 (3) + ft(4) > a (l) + fl (2) + a (3) + a (4) > ai + a2 > ^ + / g 2 . 

This implies > /^(l — c). Then it must be that &( 4 ) = /?3C, and /J3C > /^(l — c). 

Now what remains is to determine the order between and b^\ We consider b^ 
first. Nielsen's theorem yields b^ + b^ < + a 1 - 8 ). By definition, we know that 
a (7) + a (8) < a3 ( X _ c ) + a4 (i _ c ). Therefore, 

&(7) + 5(8) < a3 ( X _ c ) + a4 (! _ c ) = ( Q3 + a4 )(! _ c ) < ^ 3 + _ c ) 5 

the last inequality is due to (J2J). Since 6(8) = /3 4 (1 - c), it follows that b^ < - c). 
Furthermore, we obtain &( 7 ) = /? 4 c, &( 6 ) = /3 3 (1 - c), and /? 3 (1 - c) > /3 4 c. 

Finally, only /?2 (1 — c) leaves, so b^ = /^(l — c). Combining the above arguments, we 
finish the proof of inequality (jJJJ) . 

Clearly, inequality © implies that 

A < c< ( A A 1 (7 ) 



ft + ft I ft + ft ' ft + ft 

This is needed in the remainder of the proof. 

Remembering the order of B has been found out, it enables us to calculate easily 
Si=i for each I. The only rest problem is how to calculate Xa=i a ^ ■ To ^ ms en< ^' we 
need the following simple lemma: 

Lemma 2.1 Assume A = {a±, . . . , a n }, B = {b\, . . . , b n }. Sort B in decreasing order and 
denote the resulted sequence by 6« > 6( 2 ) > • • • > &( n ) . Then A ~< B if and only if for 
1 < I < n, 

l 

max V ai < V &« (8) 

A'CA,|A'|=Z ^ z -' 
a;gA' i=l 

with equality when I = n. 

Proof of Lemma: The "if part is obvious. For the "only if " part, we sort A in decreasing 
order and denote the resulted sequence by aW > > ■■■ > a^ n \ Then A -< S if and 
only if for 1 < z* < n, 

i=i j=i 

l 

It is easy to see that = max ai, so the lemma holds. 

^ A'CA,\A>\=1 ^ 

1 = 1 ~ CLi^A' 
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Proof of Theorem 2.1 (continued): Now the above lemma guarantees a quite easy way to 
deal with Yli=i enumerating simply all the possible cases. For example, + = 
a\c + — c) or a\c + 02c, i.e. + = max{«ic + «i(l — c), a\c + 012c}. The 
treatments for Ya=i a ^ > • • • > Yl^=i a ® are the same. What we still need to do now is to 
solve systematically the inequalities of Yli=i a — Ei=i ^ (1 — ^ — 8). We put this 
daunting but routine part in the Appendix. □ 

The above theorem presents a necessary and sufficient condition when a 2 x 2 catalyst 
exists for a transformation from one 4x4 state to another. Moreover, it is also worth 
noting that the theorem is indeed constructive. The second part of it gives all 2 x 2 
catalysts (if any) for such a transformation. To illustrate the utility of the above theorem, 
let us see some simple examples. 

Example 2.1 This example is exactly the original example that Jonathan and Plenio [0] 
used to demonstrate entanglement catalysis. Let \ip\) = (0.4,0.4,0.1,0.1) and \ip2) = 
(0.5,0.25,0.25,0). Then 

max i a 1 a ' 1 ~ H f = max { - 6 ' 1 - 2 /3) = 0.6, 

I 02 + 03 03 ~ «3 J 

min ( 131 , l3l ~ a n \ 1 = min{5/8, 2/3, 1 - 0} = 0.625. 

[ Oil + 02 «2 — P2 «3 + «4 J 

Since 0.6 < 0.625, Theorem 2.1 gives us a continuous spectrum \<f>) = (c, 1 — c) of catalysts 
for and \ip2), where c ranges over the interval [0.6,0.625]. Especially, when choosing 
c = 0.6, we get the catalyst \(p) = (0.6,0.4), which is the one given in j^j. 

Example 2.2 We also consider the example in jTO]. Let \ipi) = (0.4,0.36,0.14,0.1) and 
|^ 2 ) = (0.5,0.25,0.25,0). The catalyst for |^i) and \ip 2 ) given there is 4> = (0.65,0.35). 
Note that 

(ai+a 2 -0t a A -0 A \ -\nl-\-W n ko 

max < , 1 — — > = maxj0.52, 1 — 10/111 = 0.52, 

I 02 + 03 03 ~a 3 j 

min ( 131 / 1 ~ a 3 \ l ^—X = min{25/38, 10/11,1 - 0} = 25/38, 



«1 + «2 «2 — 02 Ct3 + "4 . 

and 0.52 < 0.65 < 25/38, Theorem 2.1 guarantees that \4>) is really a catalyst; and it 
allows us to find much more catalysts \<p) = (c, 1 — c) with c G [0.52, 25/38]. 



3 An efficient algorithm for deciding existence of catalysts 

In the last section, we was able to give a necessary and sufficient condition under which a 
2x2 catalyst exists for an transformation between 4x4 states. The key idea enabling us 
to obtain such a condition is that the order among the Schmidt coefficients of the tensor 
product of the catalyst and the target state in the transformation is uniquely determined 
by Nielsen's Theorem. However, the same idea does not work when we deal with higher 
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dimensional states, and it seems very hard to find an analytical condition for existence 
of catalyst in the case of higher dimension. On the other hand, existence of catalysts 
is a dominant problem in exploiting the power of entanglement catalysis in quantum 
information processing. Such a dilemma forces us to explore alternatively the possibility 
of finding an efficient algorithm for deciding existence of catalysts. The main purpose is 
to give a polynomial time algorithm to decide whether there is a k x k catalyst for two 
incomparable n x n states |Y>2), where k > 2 is a fixed natural number. 

To explain the intuition behind our algorithm more clearly, we first cope with the case 
of k = 2. Assume = (cci, . . . , a n ), and \ip2) = ■ ■ ■ , fin) are t wo n x n states, and 
assume that the potential catalyst for them is a 2 x 2 state (p = (x, 1 — x). The Schmidt 
coefficients of \ipi)\4>) and IV^)!^) are then given as 

A x = {atix, a 2 x, ... , a n x; ai(l-x),..., a n (l - x)} 

and 

B x = {f3 1 x,(3 2 x, . . . ,p n x;Pi(l -x),.. . ,(3 n (l - x)}, 

respectively. Sort them in decreasing order and denote the resulting sequences byaW(x) > 
a^{x) >•••> a( 2n \x) and b^(x) > U 2 \x) >■■■> b^ 2n \x). By Nielsen's theorem we 
know that a necessary and sufficient condition for \ipi)\(j)} — > |"02)|<^) is 

i i 
^aW(x)<^6«(x) (/ = l,...,2n). 

i=i i=i 

Now the difficulty arises from the fact that we do not know the exact order of elements 
in A and B. Let us now consider this problem in a different way. If we fix x to some 
constant Xq, we can calculate the elements in A, B and sort them. Then if we moves x 
slightly from xq to xo + e, the order of the elements in A (or B) does not change, except 
the case that x goes through a point x* with Qj(l — x*) = ctjX* (or /3j(l — x*) = (3jX*), i.e. 
x* = - a ^ a : (resp. x* = p.+p. ) for some i < j. This observation leads us to the following 
algorithm: 
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Algorithm 1 


1. 




2. 


Sort {pij}U{Jij} in nondecreasing order, the resulted sequence is denoted by 7W < 

7 (2) <.'..< 7 (" 2 -«) 


3. 


7 (0) <_ Q_5^(n 2 -n+l) ^_ j 


4. 


For j = to n 2 — n do 


5. 


-yW-pyti+i) 

C <— 2 


fi 


TlpfpnYi itip tnp nvn pt dt p1 pin pntc! m ,4 qtiH yp^'npp'f'i vpIv 

-L/ C uCl 111111C LllO VJ1 1_1C1 *J1 OlOlllOlluQ 111 Ou\.WX -*-*C 5 ICijUCvII V ^-'-'v 


7. 


Solve the system of inequalities: 

f ELi^O^EU^) (/ = l,..-,2n) 

\ 7M < x < 7( ,+1 ) 


8. 


OUTPUT: Catalysts do not exist, if for all i £ {0, 1, . . . , n 2 — n}, the solution set 
of the above inequalities is empty; catalyst (x, 1 — x), if for some i the inequalities 
has solution. 



It is easy to see that this algorithm runs in 0(n 3 ) time. In jl(Jj . an algorithm for the 
same purpose was also given, but it runs in 0(n\) time. 

By generalizing the idea explained above to the case of k x k catalyst, we obtain: 



Theorem 3.1 For any two n x n states = (a±, . . . ,a n ) and \ip 2 ) = (Pi, ■ ■ ■ ,Pn)> the 
problem whether there exists a k x k catalyst \<p) = (x%, . . . for them can be decided in 
polynomial time about n. Further more, if there exists a kx k catalyst, our algorithm can 
find all the catalysts in 0(n 2fc+3 ' 5 ) time. 

Proof. The algorithm is similar to the one for the case k = 2. Now the Schmidt coefficients 
of 1^1)10) and \ip2)\((>) are 

A x = {aixi, . . . ,a n xi;aix 2 , . . . ,a n x 2 ; . . . ,a n x k } 

and 

B x = {/3lX!, f3 n XD 01X2, ■ ■ ■ , (3 n X2; • • • , PnXk}- 

If we move x in the k— dimensional space the order of the elements in A x (or B x ) will 
change if and only if x goes through a hyperplane a^x^ = ctj 1 Xj 2 (fl^x^ = flj^Xj^) for some 
i\ < j\ and i 2 > j 2 . (Indeed, the area that x ranges over should be (k — 1)— dimensional 
because we have a constrain of Ei=i x i = !■) So ^ rs ^ we can write down all the equations 
of these hyperplanes 

T = {a h x i2 = a h x j2 \ii < ji,i 2 > j 2 } U {(3 h x i2 = (3j x x j2 \ii < ji,i 2 > j 2 }, 
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where |T| = 2(g)(2) = 0{n 2 ). In the k— dimensional space M fc , these 0(n 2 ) hyperplanes 
can at most divide the whole space into 0(0{n 2 ) k ) = 0{n 2k ) different parts. Note the 
number of parts generated by these hyperplanes is a polynomial of n. Now we enumerate 
all these possible parts. In each part, for different x, the elements in A x (or B x ) has the 
same order. Then we can solve the inequalities 

l 1 

i=l i=l 

and check the order constrains by linear programming. Following the well-known result 
that linear programming is solvable in 0(n 3 ' 5 ) time, our algorithm runs in 0(n 2fc+3,5 ) 
time, it is a polynomial time of n whenever k is a given constant. □ 
Indeed, Theorem 3.1 is constructive too, and its proof gives an algorithm which is able 
not only to decide whether a catalyst of a given dimension exists but also to find all such 
catalysts when they do exist. The algorithm before this theorem is just a more explicit 
presentation of the proof for the case of k = 2. 

4 Conclusion and discussion 

In this paper, we investigate the problem concerning existence of catalysts for entanglement 
transformations. It is solved for the simplest case in an analytical way. We give a necessary 
and sufficient condition for the existence of a 2 x 2 catalyst for a pair of two incomparable 
4x4 states. For the general case (kx k catalysts for n x n states), although we fail to give 
an analytical condition, an efficient polynomial time algorithm is found when k is treated 
as a constant. However, if k is a variable, ranging over all positive integers, the problem 
of determining the existence of catalysts still remains open. We believe it is NP-hard, 
since the set A x = {a\X\, . . . , a n x\; a±X2, ■ ■ ■ , a n X2] ■ ■ ■ , ce n Xk} in the proof of Theorem 
3.1 potentially has exponential kind of different orders. 

Acknowledgements: The authors are very grateful to the anonymous referees for 
their invaluable comments and suggestions that helped to improve the presentation in this 
paper. 
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5 Appendix: Proof of Theorem 12.11 

Proof of Theorem 2.1 (remaining part): We need to solve the system of inequalities 
a$ < 5ZL_i (1 < Z < 8). This is carried out by the following items: 

(1) First, we have: 

(i) < 5(1) aiC < p lC ax < (9) 

(2) The inequality + < &W + b^ may be rewritten as 

maxjaic + «i(l — c), ot\c + 02c} < 0\c + j3\{l — c) <J=^ (10) 

C < ^ , «!</?!• (11) 

ai + 02 

(3) We now consider aW + + < &W + b^ + b^ 3 \ It is equivalent to 
maxjaic + «i(l — c) + «2C, aic + 02c + 03c} < [5\c + j3\{l — c) + /?2C <J=^ 

c < < h , ft~ ai ](l2) 

a\ + a 2 + «3 - 02 ' "2 - 02 V 
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(4) It holds that 

(1) + a (2) + a (3) + a (4) < b (l) + b (2) + 6 (3) + 6 (4) 

maxjaic + «i(l — c) + a^c + 02(1 — c), a*ic + «2C + 0:3c + «i(l — c), 

o-ic + a 2 c + a 3 c + a 4 c} < 0\c + 0\{\ - c) + /? 2 c + /3 3 c 

ai + a 2 - 01 . . f A Pi - oil 1 

< c < ' 



02 + 03 ~ ~ 1 1 - 02 ~ 03 ' a 2 + a 3 - fo - 3 
(5) 

a (U + a (2) + a (3) + a (4) + a (5) < 6 (1) + b (2) + ft (3) + ft (4) + & (5) 

a (6) + a (7) + Q (8) > fe (6) + 6 (7) + ^(8) 

min{a2(l — c) + 03(1 — c) + 04(1 — c), a 3 (l — c) + 04c + 04(1 — c)} 

>/? 3 (l-c) + /3 4 c + /3 4 (l-c) 

a 4 - /3 4 /3 4 
1 — < c < 1 



/3 3 - a 3 a 2 + «3 + "4 - 03 

(6) 

Eli < Eti & w 

a (7) + Q (8) > b (7) + b (S) 

min{a 3 (l — c) + 04(1 — c), a 4 c + a 4 (l — c)} > /?4C + /?4(1 — c) 

At 
03+04 ' 

(7) We have 



i=l 1=1 

Combining Eq. (0 19I16|) we obtain 

01 03 01 01 01 -CXI 01 



c < 



01 + 02 03+ 04 OC1 + U2 OL\ + a 2 + a 3 - 02 «2 - 02 1 -02 - 03 



-i 1 £2 1 



«2 + «3 — /?2 — /?3 0L2 + a 3 + «4 — /? 3 a 3 + a 4 

and 

. ,ol\ + oli-0\ a 4 - 04 
c > < . , 1 



02 + 03 03 ~ "3 . 

Since 

01 > Oil > a 2 > 02 > 03 > «3 > «4 > 04, "1 + «2 > 01 + 02, 



if 02 + 03 — 02 — 03 < 0, this term is useless. 
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it follows that 

01 > 01 



0i + 02 atx + «2 

03 1 _ Pi 1 04 



03 +04 03+ 04 « 3 + «4 ' 

01 < 01 



ai + a>2 «i + «2 + ("3 — 02) ' 

01 c 01 01 01 



ai + a 2 01 +02 01 + 04 1-02- 

1 04 > j ft 



'3 



and 



«2 + «3 + "4 — A "3 + «4 



01 -a-t > 01 ~ &i 



Oi2 + «3 - 02 ~ 03 a 2 ~ 02 

This indicates that there are six useless terms in Eq. (jl7|) . so we can omit them. Now we 

get 

' ai + OS? - 01 a 4 - 04 \ . f 01 01- Oil 04 

max < , 1 — > < c < mm < , — , 1 



02 +03 03 — 013 J [ai + a 2 «2 —02 «3 + OL 4 

Therefore, Eq. (JIJ is a necessary condition for the existence of catalyst. 

On the other hand, we claim that Eq. Q and Eq. @ are the sufficient conditions. 
Indeed, if we choose a c satisfies Eq. (|2,1|) . then c satisfies Eq. (|17|) and ((T8|) . From 
Eq. CTTKl) we know that Yh=i « W ^ EEi^®, i-e. |V»1>|0> -» l^2>|0> under LOCC. This 
completes the proof. □ 
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